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J-H ^ Sakurai's Object (V4334 Sgr), first detected in 1996, is the cen- 
. 5r , tral star of a planetary nebula (CSPN), that experienced a Very The analysis of the very rapid evolution of the spectral en- 
Late Thermal Pulse (VLTP), a helium-shell flash while on the ergy distribution (SED) of Sakurai's Object (Eyres et al. 1999; 
white dwarf coolin g track that ca n influence the late evolution Kerber et al. 1999; Tyne et 2002; Kaufl et al. 2003; Eyres et al. 
of low-mass stars ( [Herwig 200 1| ILawlor & MacDonald 20031 2004; Evans et al. 2006; Worters et al. 2008) has suggested that 
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ABSTRACT 

Context. In 1996, Sakurai's object (V4334 Sgr) suddenly brightened in the centre of a faint Planetary Nebula (PN). This very rare 
event was interpreted as the reignition of a hot white dwarf that caused a rapid evolution back to the cool giant phase. From 1998 on, a 
copious amount of dust has formed continuously, screening out the star which has remained embedded in this expanding high optical 
depth envelope. 

Aims. The new observations, reported here, are used to study the morphology of the circumstellar dust in order to investigate the 
hypothesis that Sakurai's Object is surrounded by a thick spherical envelope of dust. 

Methods. We have obtained unprecedented, high-angular resolution spectro-interferometric observations, taken with the mid-IR in- 
terferometer MIDI/VLTI, which resolve the dust envelope of Sakurai's object. 

Results. We report the discovery of a unexpectedly compact (30 x 40 milliarcsec, 105 x 140 AU assuming a distance of 3.5 kpc), 
highly inclined, dust disk. We used Monte Carlo radiative-transfer simulations of a stratified disk to constrain its geometric and phys- 
ical parameters, although such a model is only a rough approximation of the rapidly evolving dust structure. Even though the fits are 
not fully satisfactory, some useful and robust constraints can be inferred. The disk inclination is estimated to be 75° +3° with a large 
scale height of 47+7 AU. The dust mass of the disk is estimated to be 6 x 10~^Mo. The major axis of the disk (132°+3°) is aligned 
with an asymmetry seen in the old PN that was re-investigated as part of this study. This implies that the mechanism responsible for 
shaping the dust envelope surrounding Sakurai's object was already at work when the old PN formed. 

Key words. Planetary Nebulae-Individual: Sakurai's object (V4334 Sgr) Techniques: interferometric; Techniques: high angular 
resolution;Stars: AGB and post- AGB; Stars: circumstellar matter; Stars: mass-loss 

1 ■ Introduction 1917 and appears today to be still embedded in a disk-like dust 

structure ( [Clayton et al. 2006|IHinkle et al. 20081) . 



[Hajduk et al. 2005 1). The extended very faint planetary neb- the average dust grain size has been increasing and probably the 

ula (PN) bears witness to the previous evolution of this star, rate of dust formation has been increasing as well. The results 

confirming that the latest large mass ejection during the PN were mainly based on modeling of radiative transfer in a spher- 

phase occurred several thousands years ago (Kerber et al. 1999). ical dust shell. Tyne et al. (2002) suggested that the source was 

This new 'final flash', which began 12 years ago, returned expanding and might be resolvable by 8- 10m class telescopes in 

the star very briefly to the Asymptotic Giant Branch (AGB) the mid-IR. Kerber etal. (2002) discovered of a fast outflow from 

stage, explaining why these sources are often called 'Born- the source and suggested a bipolar morphology for the fast mov- 

again' objects. As such, a final flash is astronomically very ing gas seen in optical spectra. Recently, it has been suspected 

brief (only a few tens of years), so that the observation of that there is a strong asymmetry in the circumsteUar sheU around 

such an event is rare. Only the event experi enced by V605 Aql Sakurai's Object, but t he dust sheU itself had yet to be resolved in 

can be considered as direc tly comparable ( jClayton et al. 2006[ the optical or near-IR (lEvans et al. 20061 [vanHoofet al.(2007)l ). 



Clayton & De Marco 1997| ). V605 Aql underwent a final flash in 



Observations using the Very Large Telescope Interferometer 
Send offprint requests to: O. Chesneau are presented in Section 2. In Section 3 we derive some physical 

* Based on observations made with the Very Large Telescope parameters of the dust shell using 2D radiative-transfer models 
Interferometer at Paranal Observatory under program 079.D-0415 and we discuss the results in Section 4. 
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Table 1. Observing log 



Table 2. Model parameters. 



OB 


Time 


Base 


Projected baseline 










Length 


PA 










[metre] 


[degrees] 


Sak-1 


2007-06-29T01 


U2 


-U3 


41.9 


14.6 


Sak-2 


2007-06-29T02 


U2 


-U3 


45.9 


37.1 


Sak-3 


2007-06-29T06 


U2 


-U3 


44.4 


48.4 


Sak-4 


2007-06-30T01 


U3 


-U4 


51.3 


101.0 


Sak-5 


2007-06-30T06 


U3 


-U4 


54.7 


124.3 


Sak-6 


2007-06-30T08 


U3 


-U4 


39.2 


159.5 



Calibrators: HD152334 K4III 3.99±0.07mas, HD163376 MOIII 3.79±0.12mas, 
HD169916 Kim 3.75±0.04mas, HD177716 KlIII 3.72±0.07mas. 



2, Observations 

The source was observed in June 2007 with IMIDI (Leinert et 
al. 2003, Ratzka et al. 2003) the mid-IR recombiner of the Very 
Large Telescope (VLT). The IMIDI/VLT Interferometer operates 
like a classical IMichelson interferometer combining the mid- 
IR light (N band, 7.5-13.5 jim) from two VLT Unit Telescopes 
(8.2m, UTs). A typical IMIDI observing sequence was followed, 
as described in Ratzka et al. (2007). IMIDI provided single-dish 
acquisition images with a spatial resolution of about 250 mas 
at 8.7 yum, a flux-calibrated spectrum at low spectral resolution 
(R=30) and 6 visibility spectra from the source. The observa- 
tions were performed in High_Sens mode, implying that the pho- 
tometry of the source is recorded subsequently to the fringes. 
The visibility errors range from 8% to 15%. The accuracy of the 
absolute flux calibration is better than 10%. The log of the ob- 
servations is given in Table [T] We used two difl'erent IMIDI data 
reduction packages: IMIA, developed at the IMax-Flanck-Institut 
fur Astronomic, and EWS, developed at the Leiden Observatory 
(]V[IA-hEWSEl,ver. 1.5.1). 

The IMIDI images show that the core is unresolved at 8.7 
yum, implying that its angular diameter must be less than 150- 
200 mas. The mean level of the visibility spectra is low, indi- 
cating that the source is well resolved by IMIDI/VLTI (typical 
scale 30 X 40 mas), but the spectral modulation indicates a com- 
plex source (see Fig|2l). Significant diff'erential phases (±10-20° 
peak-to-peak) were observed for all baselines except Sak.4 and 
Sak.5 (Table[T]), which exhibit phases close to zero (±5°). 

The data are complemented by two Spitzer/IRS spectra ob- 
tained in April 2005 and IMay 2007 (Ev ans et aL 2006, Evans et 
al. 2009, in preparation), and a continuum- subtracted [O III] im- 
age of the surrounding FN acquired in October 2002 with the 
instrument VLT/FORSl (Hajduk et aLlOOS ). Fig H] shows the 
Spitzer and IMIDI spectra which are characterized by a steeply 
rising flux toward the long wavelengths. The [OIII] image of the 
FN is shown in Fig [3l The appearance of the nebula is almost 
circular, with weak but significant spatial structure. 

3. Physical parameters of the disk 

The large, spectrally dependent variations of the IMIDI visibili- 
ties could be interpreted in the context of double point or ring 
structures. Simple models fail to account for the complex visi- 
bilities observed, although the diff'erence between the mean level 
of the visibilities from the baseline UT2-UT3, compared to those 
obtained with UT3-UT4 suggest a flattened structure with an ori- 
entation of the major axis of approximately 110-140°. The dif- 
ferential phases point to a significant asymmetry in the direction 
of the minor axis, i.e. ~45°. 





Parameters 


Parameter range 


T^eff (K) 


12000 


- 


Luminosity (L©) 


10000 


- 


Distance (kpc) 


3.5 


- 


Inclination (°) 


75 


3 


PA angle (°) 


132 


3 


Inner radius (AU) 


65 


10 


Outer radius (AU) 


500 


No constraint 


a 


2.0 


0.1 


P 


1.0 


0.1 


hiooAu (AU) 


47 


7 


Dust mass (M©) 


6 X 10"^ 


3 X 10"^ 




Fig. 1. IMIDI spectrum (June 2007; thick crosses) compared with 
the Spitzer spectra obtained in April 2005 (upper solid curve; 
Evans et al. 2006) and IMay 2007 (crosses; Evans et al. 2009, 
in preparation). The smooth solid line is the flux from the best 
model assuming a distance of 3.1 kpc. 

We made use o f the conti nuum IMont e Carlo radiative trans- 
fer code ]V[C3D (Wolf 2003; see also IWolf etaL 1999"), that 
solves the radiative transfer problem self-consistently. We use a 
classical model of a stratified disk of dust (Chesneau et al. 2006, 
2007 and references therein). The dust density follows a 2D law 
(both radial and vertical): 



z 



(1) 



Available at ^http://www.strw.leidenuniY.nl/~neYec/MID^ndex.html| 



where r is the radial distance in the midplane of the disk, a is the 
density parameter in the midplane, R^, is the stellar radius and 
the disk scale height, /z(r), is given by h{r) = ho {rjR^'f, where 
ho is the scale height of reference, and p is the vertical density 
parameter. 

The absence of spectral features in the IR indicates that the 
dust is dominated by amorphous carbon grains. We assume the 
standard interstellar grain size distribution (IMathis et al. 1977) 
with spherical grain radii extending from 0.005 to 1 yum. We es- 
timate the outer radius of the disk to be 500 AU consistent with 
an expansion since 1997 at a maximum velocity of 100 km s"^. 
The initial parameters (Jeff of the source and dust mass) were 
chosen following van Hoof et al. (2007). The outputs from the 
code are described in Chesneau et al. (2006, 2007). The distance, 
luminosity, and temperature of the central star are uncertain by 
a large factor (Jacoby et al. 1998). We first fixed the distance 
to be 1.5 kpc, but no satisfactory models could be found. By 
increasing the distance, the fit of the visibility curves and the 
SED steadily improved and we achieved the best fit model for a 
distance of 3.5 kpc, a value close to the distance favored by van 
Hoof et al. (2007). We note, however, that we did not explore the 
parameter space for larger distances and that the distance chosen 
is by no means a strong constraint. Our priority in this study was 
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Fig. 2. MIDI visibilities with their error bars, compared with the best model (dashed-dotted lines, see Table O. The;^^ values of the 
fits of 16.5, 7.1, 7.8, 7.6, 5.2, 4.6. The reduced for the full data set is 8.2. 



to obtain the best possible definition of the the geometry of the 
source, by fitting the visibility curves. We were not able to find 
a model with a good match to both the visibility curves and the 
SED. The best model, presented here, is the one showing the 
highest flux although it is still too low. Some models with slight 
modifications of the parameters provided better fits to the visi- 
bilities (best;^^ ~6). Given the very short timescale to form the 
disk, and the large scale height inferred, a stratified disk might 
not be the best model to account for the density profile of the 
dust. A slowly expanding torus (IPeretto et al. 2QQ7b should be 
investigated in the future. 



4. Results and discussion 

The MIDI observations provide direct evidence for a thick, 
highly-inclined disk or torus (which therefore efficiently screens 
the central source), whose parameters are described in Table 2. 
The orientation of the major axis of the disk (134+5°) coincides 
with the asymmetry detected in the old PN (130+8°). A similar 
orientation is also mentioned in Kerber et al. (2002). Such an 
alignment between the disk and PN may not be fortuitous. In the 
case of V605 Aql, th e disk and the planetary nebula A58 share 
the same major axis (IHinkle et al. 2008 ). Of importance too, is 
the low level of asymmetry detected in the old PN compared to 
what is currently observed in the disk. Is the new mass loss more 
asymmetric or is the smaller asymmetry of the PN remnant due 
to evolution and old age? Another related object worth mention- 
ing is A 30, a PN that suffered its own VLTP a few thousand 
years ago ( Borkowski et al. 19951) . The inner nebula exhibits a 
thin equatorial disk, currently being eroded by the fast wind of 
the central star, whereas its old PN is perfectly round. 

The disk is very optically thick and at 8 yum, most of the 
light emerges from a small region above the North pole (see 
Fig 131). The scale-height of the disk/torus is very large, limit- 
ing the opening angle of the polar regions, reminiscent of the 
long-lived 'wall' structure seen in the binary post-AGB HR 4049 
[Dominik et al. (20Q3)| Spherical models of the dust shell explain 
the shift of the peak of the SED toward longer wavelengths by 
suggesting that the inner boundary has receded from the star 
(IKaufl et al. 20031 Evans et al. 2006i) . If this indicated a cessa- 



tion of mass loss and dust formation, then the optical depth of 
dust would drop very quickly, assuming that the dust is moving 
radially away from the star. The estimated size of the Sakurai's 
Object disk from our observations in 2007 is small compared 
with previous estimates, and implies that the high velocities in- 
ferred from spectral lines may not be consistent with the dust 
expansion. If the expansion velocity were consistent with the 
rapid motions observed for CO and optical lines, the disk should 
already be resolvable by 8-m class telescope in the infrared 
( |Eyres et al. 2004[IKerber et al. 20021 IWorters et al."2QQ8l) . 

The fast spectral evolution of the source toward lower tem- 
peratures, as witnessed by the Spitzer and MIDI spectra, must 
be interpreted by a growing opacity of the disk, whether by ac- 
cumulation of new material or by grain growth. Only the contin- 
uous production of new dust and/or grain growth can explain the 
shift of the SED of the disk toward longer wavelengths. A more 
rapid expansion of the disk as an explanation for the SED evo- 
lution has been ruled out by the smaller disk size reported here. 
Since the dust on the inner rim of the disk is close to the sublima- 
tion temperature, and the outer radius is not well constrained by 
the N-band measurements, a dust expansion rate cannot be ac- 
curately estimated. Keplerian motion (i.e. with a very low radial 
expansion) of the dust in the disk cannot be excluded. 

The discovery of a flattened structure such as a disk is of 
great importance to the future interpretation of the outburst of 
Sakurai's Object and the VLTP phenomenon, as much of the 
previous work is based on the assumption of a spherical ejec- 
tion. Disk models may do a better job of fitting the observed 
parameters. A well-known consequence of applying a ID radia- 
tive transfer code in the dust to a 2D astronomical source is that 
the dust is put much farther than in reality, biasing strongly all 
the inferred parameters. 

The formation of a disk/torus could be explained by two 
classes of models: those involving an intrinsic asymmetry of the 
ejection, and those involving a binary companion. The models 
must also account for the asymmetry seen in the surrounding PN 
which implies that the same asymmetry has been acting on the 
mass loss for 10000 yr. Binarity can naturally explain the disk- 
like geometry and also the fast wind features (such as the He I 
/110830 line and the hole at P.A.~45°). The companion would 
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Fig. 4. Flux distribution of the model at 8, 10.5 and 13 /im. A small amount of light emerges at 8 yum while the flux comes from 
deeper in the disk at 13 jim. 





Fig. 3. Continuum subtracted [OIII] contour plot of the old PN 
recorded in 2002. The region oriented at P.A.=130° is brighter 
and a hole is visible in the North-East part of the roughly spheri- 
cal nebula, with an half-aperture angle of -45° . North is up and 
East is left. 



likely be a very low-mass object which has survived the previ- 
ous evolution of the system (Mulla lly et al. 2007 ^. The fact that 
the asymmetry of the disk is much more pronounced in the re- 
cent VLTP event could be a consequence of the smaller mass of 
the ejecta. It is important to stress that in this scenario the origi- 
nal ejection might have been spherical, but with the asymmetry 
being enhanced by the influence of the companion. The short 
timescale for the disk formation is a severe constraint and may 
imply that the orbital period of the companion is short (i.e. typi- 
cally a few months). Given their similarities, any solution to the 
asymmetry of the Sakurai's Object dust envelope will also ap- 
ply to its twin, V605 Aql ( Clayton et al. 2006). In both cases, we 
seem to have a thick disk of dust preventing the star from being 
seen directly. In both cases too, some structure in the surround- 
ing PNe are aligned with the axis of the disk. 

The small rate of expansion of the disk is a strong ar- 
gument for a binary interaction. A VLTP happens in 10-20% 
of post-AGB stars. A companion, originally far from the pri- 



mary (-5-10 AU) so to avoid mass transfer and a common 
envelope phase could be responsible for the outflow shaping. 
During the initial AGB phase, the spherical wind from the 
primary star interacts with the companion and is focused to- 
ward the equatorial plane. Initially the asymmetry is strong. 
However, after 10"^ years, the asymmetry slowly disappears 
as the nebula expands into the ISM. When the VLTP oc- 
curs, the resulting mass-loss also interacts with the compan- 
ion and is once again focused towards the equatorial plane. 
Given the mass ejected, the companion might also migrate in- 
ward and may provide enough angular momentum to stabilise 
the ejecta in the form of a t orus (Nordhaus & Blackman 2 0061 
IPeretto et al. 2007[|Edgar et al. 2008) . Given the disk inclination 
and high optical depth, the detection of such a companion will 
be a difl&cult challenge for the foreseeable future. 
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